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ABSTRACT 



A control system based on sliding mode control and the linear quadratic 
regulator is designed to stabilize the straight line motions of towed vessels. The 
control technique is through athwartship movement of the towline attachment 
point on the towed vessel. Control design is based on the linearized sway and 
yaw equations of motion. Numerical simulations for both the linearized and the 
nonlinear system are performed and demonstrate the added robustness of the 
control technique employed. 
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I. INTRODUCTION 



The horizontal plane stability of towed vessels is a very important field. If 
the motion of a vessel under tow is unstable, a very dangerous situation exists; 
the motion of a several thousand ton vessel could endanger the ship towing it 
[Refs. 1, 2, 3]. A method needs to be developed not only to ensure the stability of 
the towed vessel's motion, but to control it [Refs. 4, 5]. If the motion of the towed 
vessel can be controlled, the next step is to optimize its control, both for safety 
and fuel considerations. Usually the only means of control of a vessel under tow 
is through the towline. Under these conditions, stability of the vessel is marginal; 
in some cases, the vessel is actually unstable. 

This thesis will develop the equations of motion for a vessel in tow. It will 
also state and describe the unique forces associated with a vessel being towed. 
Finally, it will describe how state-space theory is being applied to control the 
towed vessel's motion through the movement of the attachment point of the 
towline along the beam. 

A control system based on sliding mode control and the linear quadratic 
regulator for the linearized equations of motion will be presented, with results in 
non-dimensional form. The same control system for the non-linearized equations 
of motion will also be presented, along with the results. Comparison of results for 
the linearized and non-linearized equations of motion will help determine the 
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robustness of the system. Matrix^ will be used to calculate the required gains for 
the system, and the simulation will be done using a Fortran program. 
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II. EQUATIONS OF MOTION 



The two equations of motion will be developed using a body-fixed frame of 
reference. The origin of this reference frame is at the center of gravity of the 
towed vessel. In general, the vessel has six degrees of freedom: 

1. surge along the Xg-axis 

2. roll about the Xg-axis 

3. sway along the yg-axis 

4. pitching about the yg-axis 

5. heave along the Zg-axis 

6. yaw about the Zg-axis 

This thesis will only deal with motion in the horizontal plane. Also, a constant 
surge velocity along the Xg-axis will be assumed. Other assumptions include: 

1. no towed vessel /towing vessel interaction 

2. massless, inextensible towline 

3. no wind, wave, or current disturbances 

4. the body is symmetric 
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Since motion is confined to the horizontal plane, the applicable equations are the 
sway and yaw equations. These equations are standard Principles of Nayal 
Architecture, or PNA, equations. They are as follows: 

Sway: 



m[v + ur +wp +x^(pq +f) - + r^) +z^(qr -p)] 

= ( W - B ) COS0 cosb + Yf 

Yaw: 



(la) 



m[x^{v + ur +w) ~ ( m -vr + wq)] 

= ( xw - XpB ) COS0 cos<}) + ( y^ W - y^B ) sin0 +N^ 

If we make the aboye assumptions in terms of the yariables in the equations: 



1. No yertical motion; w=0 

2. W=B 

3. Small angle motion (linearized simulation only) 

4. No pitch rate; also negligible roll; 0=0, 0=0, (})=0 

5. Forward speed equals nominal speed; u=U 

6. Center of Buoyancy=Center of grayity; Xg, yg=x^ yp 
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Then the sway and yaw equations become; 



m{v + ur +x f) = Y. (2a) 

g f 

I r +mx (v + ur) = N. (2b) 

2 g f 

The hydrostatic forces on the right hand side of the equations are then expanded 
in a standard Taylor series; 



y = Yt; + yr + yu + y.f + y , 

F V T V r towline 



N^ = N V + N r + N.v + N.f + N, .. 

F V r V T iovAine 



where Y^ = ^ | , etc. 



(3a) 

(3b) 



Note that the rudder terms, Ng and Yg have been omitted, as the towed vessel's 
rudder is amidships; there is no applied force due to the rudder. The only applied 
force is due to the towline. Substituting (3a) into (2a) and (3b) into (2b) and 
rearranging the terms yields: 

Sway: 

(m -Y )v +( -Y. )f = Yv*(Y -mu )r + Y, .. (4a) 

V 'r' V T totdine 

Yaw: 



( 



1 



-N.)f -N v=N v+Nr 

r ' V V r 



touAine 



(4b) 
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which can then be rearranged into working form as 



Y V +(Y. -tn)v + {Y. - mw )\i/ + Y.vi/ = Y, . (5a) 

V 'r ' N/ ^ towlme 

Ny+Nt>+N,\ir + (N. -7 )\i/ = N, , (5b) 

V V V z T toxvline 

where r=\}/, x-=\iX. All that remains is to describe the towline forces in greater 
detail. 



A. TOWLINE FORCE 

In Figure 1, P is the attachment point of the tow. Xp and yp denote the 
distance from the towed vessel's center of gravity, again using the standard PNA 
body-fixed coordinate system. From the diagram, we can see that the horizontal 
force associated with the towline in the sway direction is: 

T sin ( Y + \|/ ) 

The moment causing the vessel to yaw is: 

-Tx^ ( sin (y + V ) ) - Yi/p cos ( Y + V ) 

Assuming y and \|/ are small angles: 



sin(Y + v)“Y'^V“ 



y + + 

L 



+ \|/ 




y 
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Path 



o 




Figure 1. Towed Vessel Coordinate Description. 



Also, 



if =v + u sin\j/ = i; + u\j/ 
v=y -u^-,v =y -MV|/ 
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Finally, substituting for v, in equations (5a) and (5b) and rearranging terms; 
Sway: 



. T 



-(Y.u - )\\f 

~[Yu*T 



1 



(6a) 



T'y 

\i/= p 



Yaw; 






L 

f 



-[Nu^Tx 

V p 



\ 






l¥ 



(6b) 



T 

= — xy +Ty 

p^p jp 



These two equations are motion describe the towed vessel's motion. Note that the 
coordinate y describes the lateral distance from the towing path; yp describes the 
distance of the towline attachment point from the centerline of the towing vessel. 
It is for yp that the control law is being developed. 
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B. ARRANGEMENT INTO STATE-SPACE 



Let X]=y, X2=v, x.=y, then, substituting into equations (6a) and 6b) 



(r -m)x3 + y^X4 



= -Lx, +[Ym +T 
L ' ^ 



X 

1 

L , 



1 ^ 2 - 



(7a) 



Nx, + (N. -I )x. 

V 3 ' y z ' 4 



= — X x+[N u + Tx 

L p ^ ^ p 






( 7 b) 



x^ - N X. + (N u -N ) x+ — X y +Ty 



From the definition of the state variables, it can be seen that 

i, =*3 (8a) 

and 

X = X (8b) 

2 4 

These comprise the first two of four state equations. The second two state 
equations must come from equations (7a) and (7b) above. Note that these 
equations are coupled in X3 and X4. Algebraic elimination yields the final two 
equations of motion. Multiplying equation (7a) by ( and equation (7b) 

by ( -Y^) yields an equation for Xj. 
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[iN.-l ){Y. -m) -N y.]L 

= 1 [ ( N. - / ) - X y. ] X, 

V 2 py’’! 





( > 






X 




Yu^T 


1 




V 


L 

^ J 


_ 



(N. -/ ) 

y 2 



-y, 



N u + Tx 



X 

1 

L , 



)x. 



+ [ N y. - y ( N. - / ) ] X, 

+ [(y.M -y, (N. -I ) -y. (N.u -N, ]x, 

“•'p y'y z ' y'p 



+ T 



(N. - / ) 

V z - y 






(8c) 



Similarly, multiplying equation (7a) by ( -N ^ ) and equation (7b) by ( - tti ) 

and adding yields an equation for . 
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[(N. -1 )(Y -m)-N.u.]x^ 

y z V pv/y -*4 



= ±[xAY^-m)-NJx^ 

J P ^ VI 

+ {( V. - m ) 



N u + Tx 

V p 



\ 

X 

1 

L 



-N. 



Yu^T 



X 

\*-L 

L 






)x. 



(8d) 



*[N Y -N (Y.-m)]x, 



V V V V 



+ [ -N. mu +N,m -N , Y. +N. v, )x. 

*■ V \l/p p-7\U'4 





f ^ 

X 


N. 


(Y.-m) 


1 


V 


V 


A 


~~L 



+ T 



Equations (8a), (8b), (8c) and (8d) can also be expressed in state-space matrix form 







0 


0 


1 


1 

0 






1 

0 

1 


^2 




0 


0 


0 


1 


^2 




0 


^3 


— 


^31 


«32 


‘^33 


«34 


^3 


+ 




^'4 






®42 


«43 


«44^ 


^4. 




K_ 



(9) 



Where, 



D = (Y,-m) -N^Y^ 






^32 







( \ 








= D'^ { 


Y^u + T 


+ 


1 

1 

H 

N 


N^,u + TXp 


+ 
i — 1 

^ 
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a, 3 =D-- [ N,y^- y, ( % - I, ) ] 

83, = D-- [ ( y,u - y^ ) ( - I J -Y^(N^u-N^) ] 

a,,=D-- ^ [ X, ( y,-m) -A?J 



a,.=D--{ ( y^-m) 













( \ 




N u Tx 

V p 


1 ^ J 






y^u + T 




i 



} 



a,, = D-- [ N,y^ - N^y, + ( - N,u ) m ] 



r 


1 

H 

N 

1 


h-» 

+ 




L L ^ 


1 ^ JJ 







f 


1 ^S, 


-r 


1 

> 

1 


1 ^ . 


J-^J 
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III. CONTROL SYSTEM DESIGN 



A. SYSTEM AUGMENTATION 

Equation (9) defines the state-space form of the equations of motion. One 
more state must be defined and added to this equation; that of the ordered, or 
commanded, yp. This variable shall be called yp^. The usual response of yp with 
respect to time should look like a first order system, as in Figure 2: 





time 
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The equation governing this response is as follows: 



TY+y=Y 

p p yp pc 

-1 1 

or V = — V — V 

J p rr> J p rp J pc 

P P 

where Tp is the non-dimensional time constant for the towline attachment point 
motion control system. If yp^ is now considered the input to the system, and yp 
is considered a state variable, the matrix equation becomes: 



1 

1 




-1 

y 


0 


0 


0 


1 

o 


yp 






p 

0 


0 


0 


1 


0 




^2 


= 


0 


0 


0 


0 


1 


^2 


^3 




h 


®31 


«32 


«33 


®34 


^3 






K 


«41 


«42 


«43 


«44 





1 

T 

t 

0 

0 

0 

0 



( 10 ) 



where aj,, ajj, a^, a^.^, a^, a^^, a^, bg and b 4 are as previously defined. In a 
sense, the actual distance of the towline attachment point yp is treated as an extra 
state to the system. The control system will be designed based on this augmented 
system of equations. The control system will be designed based on two principles 
of state-space theory; Linear Quadratic Regulator, and Sliding Mode Control. 



B. LINEAR QUADRATIC REGULATOR 

The linear quadratic regulator function arose from efforts to find an 
optimum means of control [Ref. 6]. Normally, a gain matrix for a particular 
control system is calculated to achieve specific closed-loop pole locations. With 
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the linear quadratic regulator, or LQR, a specific performance criterion J (or "cost 
function") is defined, with the only stipulations on the poles being that they be 
negative, or stable. This criterion J is defined as: 

J = j[x'it)Qit)xit)+u'{t)R(t)uit)]dt 

The matrices Q and R are weighting matrices. The matrix Q is a state 
weighting matrix, and R is a control weighting matrix. The gains for a control 
system can be calculated based on the defined matrices Q and R. For instance, if 
the elements of the matrix Q are small relative to R, the system will tolerate large 
errors in the final state with very little control effort. Conversely, if Q is large 
relative to R, very small errors in the state X will result, but with considerable 
control effort. 

C SLIDING MODE CONTROL 

The second aspect of state-space theory to be used in controlling the towed 
vessel is that of sliding mode control. Since the control law is based on a 
linearized set of equations of motion, a lot of uncertainty in the response exists. 
Also, some of the parameters may vary. A good example for the towed vessel 
system is the towline tension T, which will certainly vary with time. A control 
law needs to be developed that will take into accoimt both the uncertainties in the 
parameters, and any dynamics that have either not been modeled, or that have 
been linearized. A sliding mode control law can ensure both stability and 
robustness of the system, with the emphasis on robustness. The LQR gains will 
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ensure stability of the system. Sliding mode control is ideally suited to a system 
where the response oscillates between set values, such as the motion of a towed 
vessel; it uses a high speed "toggling" control law to drive the system onto a 
desired "sliding plane." 

Sliding mode control takes the standard state-space system 

X = [AJa; + [b]« ^12) 

and defines a sliding plane 

0(x) = SjATj + SjATj + SjATj + (13) 

and the coefficient s, is arbitrary. Equation (13 ) can be written as 

s^x = 0 

where 

S ~ [Sj,S2, ] 

Next, define the Lyapunov function 

y(x) = i[e(x)p 

2 

stability is guaranteed, provided V (x) is a negative definite function. Another 
way to express this is 

y(x) = 6(x) 6(x) = -V 0(x) 
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therefore, 



6 = -Ti^ sign 0 



( 15 ) 



Since 



6(x) =s^(x),^ix) =s^(x) =s^(Ax + bu ) 



so, substituting for 0 in Equation (15); 



(Ax +bu) = -if sign (Q) 



and solving for u: 



u = -(s^by^ s^Ax -T\^is^by^ sign(6) 



( 16 ) 



Equation (16) is a sliding mode control law for a generic system. 



D. THE CONTROL LAW 

All that remains is to define a control law for the augmented equations of 
motion for the towed vessel. The control law is defined using sliding mode 
control and linear quadratic regulator principles. Take the augmented system: 
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11 



12 



-.1 0 0 0 0 

T 

p 

0 0 0 1 0 

0 0 0 0 1 

^3 ^31 ^32 ^33 ^34 

^4 ^41 ^42 ^43 ^44 



J 

■ 1 




' 1 ' 






Y 






p 




0 


^2 




0 








3 




0 






_ 0 _ 



’pc 



21 



22 



( 17 ) 



and split the augmented matrix A into Aj,, A,2, Aj,, and A22 as shown. Also split 
the state vector into 






Rewrite equation (17) as 



it 2 ~ ■^21^(1 ^22^2 ^2^ 

The original augmented state vector x now equals (y„ 

Next, define the sliding plane for the towed vessel as 



0(Ar) =0(y,,y2)=y, +s’'y2=0 



(18) 



and 



6 = -q^si^M(0)=y, +s’'y2 



(19) 
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but 



l/i = ^11 Vi ^ ^ “ and y, + 

so, substituting for y. and is equation (19): 

A^^y^*b^u ^sHA^^y^ ^ A^y^) = -r\^ sign{Q) 

solving for u: 



u +s’'/l 2 i)yi ^s'^A^y^] -^T\^sign(e) (20) 

This is the system control law for the towed vessel. 

The gains [s] for the sliding mode control law from minimization of the 
linear quadratic regulator cost function J. Recall that 

J=j{[xy[Q][x]^[uy[R][u])dt <21) 

in Equation (21) the state weighting matrix Q is as follows: 

'O 0 0 o' 

0 0 0 

^ 0 0 0 0 

0 0 0 0 

where q ^2 represents an error in radians from the path of the towing vessel. We 
are choosing to weigh the state variable Xj = \\f. The control weighting matrix R 
will be defined by 
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R = (a)-^ 

where a, the maximum non-dimensional distance will be weighted. With Q 
and R defined, the cost matrix function becomes: 

] = min * Ru^)dt ^22) 

or 

/ = min J ( ^22^ ^ ^ 

The smaller the value for R is, the smaller the control effort required, but a larger 
state error in q 22 will be have to be accepted as a trade-off. Conversely, the larger 
the value for R is the more control required, with very little state error in q 2 i- 
The resulting gains are then placed into the control law for the augmented 
system of equations. Equation (20), and the response of the system is obtained. 
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IV. RESULTS 



A. MATRIXX AND FORTRAN IMPLEMENTATION 

MatrixX is an outstanding tool for understanding the response of a control 
system. A Fortran program was written to use with MatrixX. Since the coefficients 
of the [A] and [B] matrices of the state-space equation are constant with time, 
they can be put directly into the MatrixX software, and the resulting Linear 
Quadratic Regulator gains are easily calculated. Without the use of MatrixX, 
solving the linear quadratic integral is a very difficult problem in numerical 
analysis. The actual simulation of the system response was accomplished with a 
Fortran program on the VAX computer. The LQR gains obtained from MatrixX 
are put into a data file, along with the non-dimensional parameters of Tp, LL, Xp, 
T, and the initial conditions of Y, V|/, and the maximum distance Yp can travel 
athwartships. The simulation then reads this data file, computes the sliding mode 
control and plots the response of the parameters Y, Yp, \y, and Yp<. with respect to 
time. 

B. INITIAL CONDITIONS AND CONSTANT PARAMETERS 

The results are presented graphically, and are in non-dimensional form, with 
the exception of the yaw angle \j/, which is in degrees. The non-dimensionalization 
is standard Principles of Naval Architecture: the relationships describing the 
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nondimensional parameters are included in the Appendix. Two vessels were 
studied: 

1. 528 foot mariner 

2. 1066.3 foot tanker 

A summary of the hull particulars and hydrodynamic derivatives for both these 
vessels in also included in the Appendix. Four knots was used as the nominal 
forward velocity. The mariner is stable at four knots; the tanker is unstable. The 
towline tension for both vessels is taken from resistance curves at four knots, and 
then non-dimensionalized; the value used is 0.001 Xp, the longitudinal distance 
from the towed vessel center of gravity to the towline attachment point, is 
assumed to be constant. A value of 0.5 is used. The time required for Yp to 
"match" Yp^ as in Figure 2, Tp is set as 0.5. The maximum non-dimensional 
distance Yp can travel from port to starboard extremes is taken to be 0.1; therefore 
the maximum distance from an extreme to centerline is 0.05. The initial conditions 
for each data run are Y=0.5, V|/=0, and Yp=0.05 as in Figure 3. Also, the weighting 
matrices were set at q22=5° (or .087 radians) for Q, the state error weighting 
matrix, and a=0.015 for R, the control effect weighting for Yp^. Four plots were 
generated for each simulation. The y-distance from the path (Y), the towline 
attachment point lateral offset (Yp), the towline attachment point offset command 
(Yp^) and the yaw angle (y) are all plotted vs nondimensional time. Variance of 
parameters of interest are shown on subsequent runs. 
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Figure 3. Initial Conditions. 
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1. Figures 4 and 5: Standard Mariner and Tanker 

Figure 4 shows the response of the mariner with the "standard" 
parameters and initial conditions stated above; Figure 5 shows the "standard" 
tanker. Overshoot in Y and V|/ exists for both cases, although it is higher in the 
mariner. The Yp^ graph shows that significantly more effort is required to control 
the tanker, and the oscillation in Yp is much greater than for the mariner. This 
makes sense, as the tanker is inherently more unstable. Both vessels converge to 
a straightline path in approximately 60 non-dimensional seconds. As the tanker 
is more unstable, it will be used to demonstrate all other variations in parameters. 
Henceforth all references to "normal" or "standard" tanker response will be to the 
response shown in Figure 5. 

2. Figure 6: Taidcer with Y, Y Not Observable 

Figure 6 shows the vessel response if the state variables Y and Y are not 

observable and are not used in the control law. Both Y and Y converge but the 
time required is significant; over 100 seconds. In this case there is no overshoot 
in Y, but considerable overshoot in \j/. \j/, Yp^ and Yp all oscillate considerably; but 
in the end, the system does converge. 

3. Figures 7, 8: Observable Tanker, Variance of 0(x) 

Recall from Equation (13) that the sliding plane 0(x) is in part arbitrarily 
defined, as the gain Tj^ is chosen by the designer. For the standard runs, was 
set to be 1.0. Figure 7 shows the tanker response if T|^ in 0(x) is set to 0.5, with all 
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Figure 4. Mariner Response with Standard Parameters. 
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Figure 5. Tanker Response with Standard Parameters. 
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Figure 6. Tanker with Y, Y Not Observable. 
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state variables observable. The system still converges in adequate time; the 
parameters Y, Yp^, and Yp behave as in the standard tanker. There is 
approximately twdce as much overshoot in \j/, however. In Figure 8, ri^ is set to 
2.0; this brought the overshoot in \\f down to the level of the standard observable 
tanker, without altering the favorable response of the other three parameters. 

4. Figures 9, 10, 11: Non-Observable Tanker, Variance of 6(x) 

The state variables Y and Y were assumed to be not observable in 
Figures 9 and 10. In Figure 9, was set to 0.5, with disastrous results. The towed 
vessel becomes unstable. All four parameters diverge with time. In Figure 10, 
however, T|^ in 6(x) is set to 2.0, and the system response stabilizes. The overshoot 
in y is approximately as the same as the standard tanker, but the overshoot in Yp 
and Yp^ is twice that of the normal tanker. The settling time in Y is double that 
of the normal tanker, but the settling time in, Yp^ and Yp is cut to forty seconds. 
Obviously there is a value between rf = 0.5 and = 1.0 in 6(x) where the 
response turns unstable for the non-observable case. In subsequent simulations, 
is set to be 2.0. Figure 11 summarizes the response of Y and Yp for the 
observable and non-observable, with q^ = 2. The top two graphs are for the 
observable case, and the bottom two describe the non-observable case. 

5. Figure 13; Variance of the Sliding Plane Switching 

Since sliding mode is a high-speed switching, or toggling method or 
control, the speed at which the switching is to be accomplished must be defined. 
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Figure 7. Observable Tanker, 'n^=0.5. 
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Figure 8. Observable Tankers, tj^=2.0 








Figure 9. Non-Observable Tanker, ti^=0.5. 
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Figure 10. Non-Observable Tanker, ti^=2.0. 
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Figure 11. Observable, Non-Observable Comparison 
of Y, Yp for ti^=2.0. 
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This switching speed can be defined in terms of an angle. Note in Figure 12 that 
as the angle 4) approaches 90 degrees, the system is being asked to change 
instantaneously. Another way of saying this is that the system is ordering the 
towline attachment point to move instantaneously from port to starboard. ())=90 
degrees is impossible; the designer has to allow the system time to react. 

Figure 13 shows the observable and non-observable tanker response of 
Y, Yp for <|)=30, 4)=45, and <))=60. For the observable case, overshoot, the amount of 
oscillation, and settling time increase as <() increases. For the non-observable 
tanker, the response destabilizes between <])=45 and 0=60. Clearly there is a limit 
to the switching speed of the system. 0=45 was used for the standard tanker. 

6. Figure 14: Variance of LQR Weighting 

Figure 14 demonstrates the effect of varying the control effort in the 
linear quadratic regulator. As a is increased from 0.005 to 0.035, more effort is 
used to control the towed vessel motion. In Figure 14, Y and Yp responses are 
plotted for four different values of a. In all four graphs, the tanker is observable. 
Note that as a increases, the overshoot in Yp decreases, and the settling time 
remains about the same as for the standard tanker. This makes sense; more effort 
is being expended to control Yp through Yp<.. The response of Y remains essentially 
unchanged. 
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Figure 12. The Sliding Plane. 
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Figure 13. Y, Yp Variance When Changing Switching Angle 0. 
0=30°, 0=45°, 0 = 60° 
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Figure 14 . Tanker Response for Y, Yp as a Varies. 

Y Y. 

Ip/ 
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7. Figure 15: Observable Tanker, Variance of Tp 

Figure 15 demonstrates the effect of varying the time allotted for Yp to 
match for the standard tanker. The T_ used for the standard case was 0.5. 

pc p 

Varying Tp does not significantly affect the behavior unless the change is in order 
of magnitude. This result has favorable implications; if changing Tp doesn't 
radically affect the towed vessel response, a smaller motor can be used to drive 
the device moving the towline attachment point, Yp. 

8. Figure 16: Variance of Towline Length 

Figure 16 illustrates how changing the towline length will change the 
response of the tanker, both in the observable and non-observable case. For this 
simulation, the non-dimensional towline length, LL, was lengthened and 
shortened by twenty percent. If the towline is shortened, the overshoot and 
oscillation in Y and Yp will increase for both observable and non-observable cases; 
the reverse is true is the towline is lengthened. Settling time will remain the same. 

9. Figure 17: Variance of Towline Tension 

Figure 17 shows how changing the non-dimensional towline tension 
will change tanker response. If the tension in decreased, the overshoot remains 
about the same; oscillation and settling time in Y and Yp increase for both the 
observable and non-observable cases. The reverse is true as the tension is 
increased. A 20 percent change in towline does not cause a radical change in 
system response. 
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Figure 15. Observable, Non-Observable Response 
as Tp Varies. y^,, y. 
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Figure 16. Observable, Non-Observable Response 
for Variance of LL. 

LL=2.0, LL=2.5, LL=3.0. 
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Figure 17. Observable, Non-Observable Tanker: Variance of T. 
T=.0008 T=.0010 T=.0012 
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10. Non-Linear (Large Angle) Motion 

Recall from page 6 that the horizontal towline force associated in the 
sway direction is 

T sin ( y + \j/ ) 

and the moment due to the towline causing the vessel to yaw is: 

-TXp ( sin ( Y + Y ) ) - T y^cos ( y + y ) 

Previously the assumption was made that the range of motion for y and y was 
less than 30 degrees. This simulation lifts this restriction; the sine and cosine 
terms from the towline tension term are left in when running the simulation. 
Figure 18 plots the standard tanker response for both linear and non-linear 
tensions. Note that the tanker response remains stable; the sliding mode control 
law is robust enough to handle non-linearities in the towline tensions. The non- 
linear Fortran simulation program is included in the Appendix. 
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Figure 18 . Non-Linear, Linear Tankers With Standard Parameters, 



V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

This thesis explored the possibility of athwartships movement of the towline 
attachment point as a means of improving the stability of towed vessels. Newell 
demonstrated that use of a full order observer will work only when the towed 
vessel trajectory is close to that of the towing vessel [Ref. 5]. Sliding mode control, 
in conjunction with the Linear Quadratic Regulator cost function, provides a 
much more robust method of controlling the towed vessel motion. The sliding 
mode control law can ensure stability of the system, even if not all the state 
variables can be observed. The only restrictions for the numerical applications 
considered in this work are that a toggling speed such that ([) is less than an angle 
between 45 and 60 degrees is used, and that the toggling gain T|^ is greater than 
a number between 0.5 and 1.0. The LQR function allows the designer to choose 
where his emphasis is placed. By choosing whether to place the emphasis on the 
control effort, the designer can minimize the size of the motor required to drive 
the mechanical apparatus used to move the towline attachment point. If the 
designer wants to minimize steady state error, he can change the LQR function 
to reflect that emphasis. The LQR function can also be used to minimize specific 
state variables, such as in the case studied, where \|/ was chosen to be minimized. 
It should be emphasized that a "worst-case" scenario was used to demonstrate the 
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robustness of the system; a large, unwieldy tanker was used, and two of the four 
state variables were declared "not observable." Even in the worst case the sliding 
mode control law performed adequately. Towline tension and length were 
increased and decreased by 20 percent, and the system remained stable. The small 
angle motion assumption was dropped, and yet the system continued to perform 
satisfactorily. 

B. RECOMMENDATIONS 

The sliding mode control law should be robust enough to cover most of the 
unmodeled dynamics and non-linearities, but two of the assumptions made at the 
beginning of the thesis must be investigated further; the massless, inextensible 
towline, and the no wind, wave, or current (disturbances) assumptions. In reality 
the towline tension does not remain constant; consequently, the surge velocity u 
of the towed vessel will also vary. A method to take these aspects of motion into 
account needs to be looked at; integral control is one possibility. Another is to 
determine a wave and towline function, and add the two functions to the right 
side of the sway and yaw equations of motion. 

This thesis has shown that it is possible to stabilize and control the motion 
of a towed vessel with a sliding mode control law that uses linear quadratic 
regulator gains. Further analysis needs to be done on the effects of disturbances 
and a time-variant towline tension on the system. 
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APPENDIX 



TABLE I. TOWED VESSEL DATA 

Vessel 



Property 


Barge 


Tanker 


Mariner 


LBP, ft 


191.56 


1066.3 


528 




-0.00136 


-0.0009 


-0.0004444 




0.170 


0.0181 


0.00888 


V 


-0.01383 


-0.0171 


-0.00912 


y' 


-0.0153 


-0.0261 


-0.01434 


-Ty 


0 


0 


0 


r; 


0.00238 


0.00365 


0.00456 




0 


0 


0 


Ni 


-0.007285 


-0.0105 


-0.0046 


1 

N 

H 


0.00188 


0.00222 


0.00115 




-0.00128 


-0.0048 


-0.00296 
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TABLE 2. NONDIMENSIONAL TERMS 



m' = m / 
u' = u / u_ 

F 

v' = V / 

u' = u / { ul - L) 
v' = V / ( ul - L) 
i' = i /{ ul - L) ( time ) 
= I, / ( /n,L^ ) 
r'=\)/=\(// { u^/ L) 
t' =V|/ =V|/ / { ul / ) 

Xg= Xg/ 1 
yl=Yg/ L 

xi ^XJ ( m^uj V ) 
X[=XJ 



= (m,uj L) 

yI^YJ 

yI = Y^ = Y^/ (m^uj 
yi = / ( m^L ) 

nI = N^/ im^uj 
nI=N^/ { m^L ) 

N'r=N^ =N^/ ( m^u^L) 

nI =w' =N^/ { ) 

T' = T/ ( m^ul / L) ( 
1' = 1 / L 

Xp = Xp / L 
yp=Yp/ L 



Tension ) 
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MATRIX, LQR PROGRAM 



Xudot=-0.0009; 

M=0.0181; 

Yvdot=-0.0171; 

Yv=-0.0261; 

Yrdot=0.0; 

Yr=0.00365; 

Ydel =0.0278; 
Nvdot=0.0; 
Nv=-0.0105; 
Nrdotlz=-.00222 



Nr=-0.0048 



Ndel=-0.0139 
inquire LL 
inquire Xp 
inquire T 
U=1.0; 

L=1.0; 

D=Nrdotlz*(Yvdot-M)-Nvdot*Yrdot; 

al3=1.0; 

a24=1.0; 

a31 =(1 /D)»(T/LL)»(Nrdotlz-Xp»Yrdot); 

a32=(l/D)»((Yv»U+T*(l+Xp/LL))*NrdotIz-Yrdot»(Nv»U+'rXp*{l+Xp/LL))); 
a33=(l /D)»(Nv*Yrdot-Yv*NrdotIz); 
a34=(l/D)»((Yvdot»U-Yr)»NrdotIz-Yrdot»(Nvdot»U-Nr)); 
a41=(l/D)»(T/LL)‘(Xp*(Yvdot-M)-Nvdot); 

a42=(l/D)*((Yvdot-M)*(Nv»U+'rXp»(l+Xp/LL))-Nvdot*(Yv*U+r^(l+Xp/LL))); 

a43=(l/D)»(Nvdot*Yv-Nv»(Yvdot-M)); 

a44=(l/D)*(Nvdot*Yr-Nr»Yvdot+(Nr-Nvdot*UrM); 

b3=(l/D)»T*(NrdotIz/LL-Yrdot*(l+X'p/LU); 

b4 =(1 /Dn^(( Yvdot-M )»( 1 +Xp/LL)-N vdot / LL); 

A=(0,0^13A-0A0,a24;a3La32,a33,a34;a41,a42,a43,a44]; 

B=[0,-0;b3;b4]; 

C=[0,0,L0]; 

D=[0]; 

Q=10 0 0 0;0 131.3316 0 0;0 0 0 0;0 0 0 0]; 



inquire alpha 
R=[(l/alpha)»»2]; 

<ElG,K>=REGULATOR(A,B,Q,R); 

INQUIRE TP 

A1=(-1/TP,0A0,0;0A0,1,0A0,0A1;B3,A3LA32,A33,A34;B4,A41,A42,A43,A44]; 

Bl=(l/TP;0,-0,-0;0]; 

S=[1,K]; 



G1=-INV(S»B1)»S»A1; 

G2=-INV(S*B1); 

ETA2=G2 

GAIN=G1 

SIGMA=S 
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LINEAR SIMULATION PROGRAM 



REAL KYP,Kl,K2,K3,K4,L,LL,Nvdot,M,Nv;VIr,Nrdot 

Xudot=-0.0009 

M =0.01 81 

Yvdot=-0.0171 

Yv=-0.0261 

Yrdot=0.0 

Yr=0.00365 

Ydel=0.0278 

Nvdot=0.0 

Nv=-0.0105 

Nrdot=-0.00222 

Nr=-0.0048 

Ndel=-0.0139 

U=1.0 

L=1.0 

OPEN(10,nLE='LINEAR.DAT',STATUS='OLD') 
OPEN (1 l,nLE='Y.DAT',STATUS='NEW') 

OPEN (12,FILE='PSI.DAT',STATUS='NEW') 

OPEN (13,nLE='YP.DAT',STATUS='NEW') 

OPEN (14,nLE='YPC.DAT',STATUS='NEW') 

READ (10,») TSIM,DELTAT,IPRNT 
READ (10,‘) KYP,K1,K2,K3,K4 
READ (10,») SYP,S1,S2,S3,S4 
READ (10,‘) GEAR 
READ (10,*) TP,LL,XP,T 
READ (10,*) ETA2,PHI 
READ (10,*) Y,PSI,SAT 

ISIM=TSIM /DELTA! 

PI=4.0*ATAN(1.0) 

PH1=PH1*PI/180.0 

PS1=PSI*PI/180.0 

SATP= SAT 

SATM=-SAT 

XI =Y 

X2=PS1 

X3=0.0 

X4=0.0 

YP=0.0 

YPC=0.0 

D=Nrdot*(Yvdot-M)-Nvdot*Yrdot 

al3=1.0 
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a24=1.0 

a31=(l/D)»(T/LLr(Nrdot-Xp»Yrdot) 

a32=(l/D)»((Yv»U+'r(l+Xp/LL))*Nrdot-Yrdot*(Nv»U+T*Xpni+Xp/LL))) 
a33=(l /D)»(Nv*Yrdot-Yv»Nrdot) 
a34=(l/D)*((Yvdot*U-YrrNrdol-Yrdor(Nvdot*U-Nr)) 
a41=(l/Dr(T/LL)»(Xp»(Yvdol-M)-Nvdot) 
a42=(l/Dr((Yvdot-Mr(Nv»U+T*Xp*(l+Xp/LD) 

& -Nvdot»(Yv*U+T*0+Xp/LL))) 
a43=(l /Dr(Nvdot»Yv-Nv*(Yvdot-M)) 
a44=(l/D)»(Nvdot»Yr-Nr»Yvdot+(Nr-Nvdot»U)»M) 
b3=(l/D)*’r(Nrdot/LL-Yrdot*(l+Xp/LL)) 
b4=a/DrT*((Yvdot-Mr(l+Xp/LL)-Nvdot/LL) 

SIMULATION BEGINS 

DO 1 I=1,ISIM 
YPDOT =-YP/TP+YPC/TP 
XIDOT = X3 
X2DOT = X4 

X3DOT = B3*YP + A3UX1 + A32*X2 + A33*X3 + A34»X4 
X4DOT = B4*YP + A4UX1 + A42‘X2 + A43*X3 + A44*X4 

YP = YP + YPDOT»DELTAT 
XI = XI + XlDOrOELTAT 
X2 = X2 + X2DOTDELTAT 
X3 = X3 + X3DOT»DELTAT 
X4 = X4 + X4DO'TDELTAT 

CONTROL LAW 

SIGMA = SYP»YP + SrXl + S2»X2 + S3‘X3 +S4*X4 
SATSGN= SIGMA/TAN(PHI) 

IF (SATSGN.GE.(+1.0)) SATSGN=+1.0 

IF (SATSGN.LE.(-l.O)) SATSGN=-1.0 

UHAT = KYP’YP + KrXl + K2*X2 + K3»X3 + K4*X4 

UBAR = ETA2»GBAR-^SATSGN 

YPC = UHAT+UBAR 

IF (YPC.GE.SATP) YPC=SATP 

IF (YPC.LE.SATM) YPC=SATM 

Y =X1 
PSI=X2 



PRINT RESULTS 



J=J+1 

IF (J.NE.IPRNT) GO TO 1 
TIME=PDELTAT 
WRITE (V) TIME,Y 
WRITE (11/) TIME,Y 
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WRITE (12,*) T1ME,PS1»180.0/P1 
WRITE (13,*) T1ME,YP 
WRITE (14,*) TIME,YPC 
J=0 

1 CONTINUE 
STOP 
END 



NON-LINEAR SIMULATION PROGRAM 



C NONLINEAR SIMULATION - TANKER 
C SURGE NOT INCLUDED 
C 

RE AL K YP,K1 ,K2,K3,K4,L,LL,NVDOT,NV,LZ,NRDOT,NR,NVRR 
C 

XUDOT=-0.0009 
M = 0.0181 
YVDOT=-0.0171 
YV =-0.0261 
YRDOT= 0.0 
YR = 0.00365 
NVDOT= 0.0 
NV =-0.0105 
IZ =0.0 
NRDOT=-0.0022 
NR =-0.0048 
YVRR =-0.045 
NVRR = 0.0061 
UTOW = 10.0 
L = 1066.3 
RHO = 1.9905 
SB = 680625.0 
P = 9.78 
Q = 1.93 

SB = SB/(0.5*RHO*L*L*UTOW»UTOW) 

C 

OPEN dO,nLE='LINEAR.DAT',STATUS='OLD') 

OPEN dl,nLE='Y.RES',STATUS='NEW') 

OPEN (12,nLE='PSI.RES',STATUS='NEW') 

OPEN (13,HLE='YP.RES',STATUS='NEW') 

OPEN d4,HLE='YPC.RES',STATUS='NEW') 

OPEN (15,FILE='NONLINEARl.RES',STATUS='NEW') 

C 

READ (10,d TSIM,DELTATTPRNT 
READ (10,d KYP,K1,K2,K3,K4 
READ (10,d SYP,S1,S2,S3,S4 
READ (10,») GBAR 
READ (10,*) TP,LL,XP,T 
READ (10,*) ETA2,PHI 
READ (10,*) Y,PSI,SAT 
C 

ISIM = TSIM/DELTAT 
PI = 4.0*ATAN(1.0) 

PHI = PHI*PI/ 180.0 
PSI = PSI*PI/ 180.0 
SATP = SAT 
SATM =-SAT 
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X = LL-XP 
U = 1.0 
V = 0.0 
R = 0.0 
YP =0.0 
YPC =0.0 
C 

DEN=(M-YVDOT)*(IZ-NRDOT)-YRDOTNVDOT 
SIMULATION BEGINS 
DO 1 I=1,ISIM 

F2=YV»V+0.5‘YVRR»V»R»R+(YR-M)»R 
F3=NV»V+0.5»NVRR»V»R*R+NR*R 

SING=( Y+XP»SIN(PSI)+ YP»COS(PSD) / LL 
GAMMA=ASIN(SING) 

YPEX)T =-YP/TP+YPC/TP 
VDOT =(F2-'TSIN(GAMMA+PSI))/(M-YVDOT) 

ROOT =(F3-'TXP»SIN(GAMMA+PSI)-’rYP»COS(GAMMA+PSI))/(IZ-NRDOT) 
YDOT =U»SIN(PSI)+V*COS(PSI) 

PSIDOT=R 

YP = YP + YPEXDT »DELTAT 
U = U + UDOT »DELTAT 

V = V + VDOT »DELTAT 
R = R + RDOT »DELTAT 
X =X +XDOT »DELTAT 

Y =Y +YDOT »DELTAT 
PSI = PSI + PSIDOT»DELTAT 

CONTROL LAW 

SIGMA = SYP*YP + SDY + S2*PSI + S3*YIX)T + S4»PSIDOT 
SATSGN= SIGMA/TAN(PHI) 

IF (SATSGN.GE.(+1.0)) SATSGN=+1.0 
IF (SATSGN.LE.(-l.O)) SATSGN=-1.0 

UHAT = KYP’YP + K1*Y + K2»PSI + K3»YDOT + K4*PSIDOT 
UBAR = ETA2*GBAR»SATSGN 
YPC =UHAT+UBAR 
IF (YPC.GE.SATP) YPC=SATP 
IF (YPC.LE.SATM) YPC=SATM 

PRINT RESULTS 



J=J+1 

IF (J.NE.IPRNT) GO TO 1 
TIME=PDELTAT 
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WRITE (*;) TIME,Y 

WRITE T1ME,Y 

WRITE (12,») TIME,PSI»180.0/PI 

WRITE (13,*) TIME,YP 

WRITE (14,*) TIME,YPC 

WRITE (15,*) TIME,Y,PS1*180.0/PI,YP,YPC 

J=0 

1 CONTINUE 
STOP 
END 
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